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1. INTRODUCTION 

Chiral Effective Field Theory (EFT) has become a standard tool for analyzing the prop- 
erties of hadronic systems at low energy where the perturbative expansion of Quantum 
Chromodynamics (QCD) in powers of the coupling constant cannot be applied, see jl!2] . 
It is based upon the approximate and spontaneously broken chiral symmetry of QCD. 
[ Starting from the most general effective Lagrangian for Goldstone bosons (pions in the 
two-flavor case of u and d quarks) and matter fields (nucleons, A, . . .) consistent with 
the symmetries of QCD, the hadronic S-matrix elements are obtained via a simultaneous 
expansion in the low external momenta and quark masses. Goldstone boson loops are 
.£h ! naturally incorporated and all corresponding ultraviolet divergences can be absorbed at 
each fixed order in the chiral expansion by the counter terms of the effective Lagrangian. 

This perturbative scheme works well in the pion-nucleon and pion-pion sectors, where 
the interaction vanishes at vanishing external momenta in the chiral limit. In the case 
of few interacting nucleons one has to deal with non-perturbative problems. Indeed, 
perturbation theory is expected to fail already at low energy due to the presence of 
the shallow few-nucleon bound states. A suitable non-perturbative approach has been 
suggested by Weinberg who showed that the strong enhancement of the few-nucleon 
scattering amplitude arises from purely nucleonic intermediate states. Weinberg suggested 
to apply EFT to the kernel of the corresponding scattering equation, which can be viewed 
as an effective nuclear potential. This idea has been explored in the last decade by many 
authors. In the following I will briefly address some of the actual topics related to the 
EFT description of few-nucleon systems. 



2. FEW NUCLEONS AT NEXT-TO-NEXT-TO-LEADING ORDER 

The procedure suggested by Weinberg has been first carried out for two nucleons by 
Ordonez and co-workers, who derived a NN potential up to next-to-next-to-leading order 
(NNLO) in the chiral expansion and performed a numerical analysis of the two-nucleon 
(2N) system 4j. The explicit energy dependence of the effective potential derived in jl] 
leads, however, to difficulties in applications to > 2N systems but can be eliminated 
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Figure 1. np differential cross section (upper row) and vector analyzing power (lower 
row) at £/i a b = 50 MeV (left panel) and E-^ = 96 MeV (right panel). Shaded band 
refers to the NNLO result, dashed line to the Nijmegen phase shift analysis (NPSA) [13J. 
References to data can be found in [T3] . 



by certain techniques. The corresponding energy-independent expressions for the NN 
potential have been first given in ref. [5] and later in refs. |6|7j using different methods. 
They served as a basis for the NNLO analyses of the two-nucleon (2N) system |8l9j . see 
also JU| for related application, as well as three- (3N) and four-nucleon (4N) systems 
[TT] including the corresponding 3N forces. 

According to chiral power counting the dominant contributions to the effective Hamil- 
tonian for few nucleons are of the order Q°/A°, where Q ~ refers to the soft scale 
(typical momenta involved in the process) and A to the hard scale (chiral symmetry break- 
ing scale, ultraviolet cut-off) [3]. These leading-order (LO) contributions turn out to be of 
the 2N type, given by one-pion exchange (OPE) and two NN contact interactions without 
derivatives. The first corrections at next-to-leading order (NLO) are of the order Q 2 /A 2 
and still of the 2N type. They result from two-pion exchange (TPE) with the leading (i.e. 
with one derivative) tcNN and mrNN vertices and 7 independent NN contact interactions 
with two derivatives. 1 The only low-energy constants (LECs) entering the expressions 
for TPE at NLO are the nucleon axial-vector coupling and the pion decay constant 
F n . Both LECs are measured rather accurately, so that the leading TPE contribution 
is parameter-free. On the contrary, the LECs accompanying the contact operators are 
unknown and have been fixed from a fit to low-angular-momentum partial waves. 

At NNLO (~ Q 3 /A 3 ) one has to take into account the subleading TPE contributions 
given by the triangle diagram with the ttttNN vertex with two derivatives or one insertion 
of M%. The corresponding LECs are denoted 01,3,4 and have been fixed in the ttN system, 
see e.g. ^H]- The numerical values of these LECs found in ^H] as well as in several other 

1 The additional two momentum-independent contact interactions with one insertion of M£ can be ab- 
sorbed by the LO contact operators. 
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Figure 2. Nd differential cross section and vector analyzing power at = 65 MeV. 
Light (dark) shaded bands refer to the NLO (NNLO) results. Data are from |14|15j . 



analyses of ttN scattering are rather large compared to what is expected on dimensional 
reasons, see ref . [H] . Similar large values of 03,4 have also been obtained recently from the 
np and pp partial wave analysis carried out by the Nijmegen group jTH]. The large values 
of the C3 5 4 can at least be partially explained by the fact that these LECs are saturated 
by the A-excitation [20] • This implies that a new and relatively small scale, namely 
rriA — tun ~ 293 MeV enters the values of these constants in EFT without explicit 
A. The large numerical values of the q's lead to the sub leading TPE contribution to 
the NN potential which shows an unphysically strong attraction already at intermediate 
distances r ~ 1 — 2 fm when standard regularization techniques to pion loop integrals (i.e. 
dimensional or infinite momentum cut-off regularization) are applied. The unphysical 
behavior of the potential shows up most notably at NNLO in parameter-free predictions 
for D-waves (F-waves) which start to strongly deviate from the data already at > 50 
MeV (Ejab > 150 MeV) after the subleading TPE contribution is taken into account jB]. 
Higher order counter terms are needed in order to reduce the strong cut-off dependence of 
the D-wave phase shifts and to correct the F-waves, which indicates the slow convergence 
of the chiral expansion. In addition, unphysical deeply bound states arise in low partial 
waves jHJ. Although such deeply bound states do not influence NN observables at low 
energy, they might lead to certain complications in another processes like e.g. Nd (Hj and 
ird [2T] scattering. 

In order to avoid the above mentioned difficulties we have constructed in ref. the 
NNLO* version of the potential without deeply bound states. To achieve that, we adopted 
values of the LECs 03,4 much smaller in magnitude than the ones obtained from ttN 
scattering. The resulting potential leads to a good description of the low-energy NN data 
as exemplified with selected observables in Fig. ^ 

Lot of progress has also been achieved for 3N and 4N systems. Scattering and bound 
state problems in such systems can be solved in a numerically exact way using Faddeev- 
Yakubovsky theory and serve as a good testing ground for chiral forces, since most of 
the unknown parameters are already fixed in the two nucleon system. As already stated 
before, no 3N forces appear at LO and NLO. Parameter-free calculations for Nd scattering 
as well as 3N and 4N bound states at that order in the chiral expansion have been presented 
in ref. [TH]. We found a good description of various scattering observables at low energy, 
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see Fig|21 The predicted 3 H and 4 He binding energies 

BE( 3 H) = -7.53 ... - 8.54 MeV for A = 500 . . . 560 MeV, 

BE( 4 He) = -23.87 ... - 29.57 MeV for A = 500 . . . 560 MeV, 

where A refers to the momentum cut-off in the 2N force, are in a similar range as the 
ones obtained using modern phenomenological NN potentials and have to be compared 
with the empirical values BE( 3 H) = -8.68 MeV and BE( 4 He) = -29.8 MeV. 2 

Chiral 3N forces start to contribute at NNLO and are given by TPE, OPE with the 
pion emitted (or absorbed) by the NN contact interaction and 3N contact interaction. 
The TPE contribution is parameter-free. The Pauli principle together with the usual 
symmetry requirements (parity invariance, rotational invariance, . . .) lead to a strong 
reduction of the number of independent terms in the remaining part of the 3N force, 
leaving just one OPE and one contact operator [TT] . The two corresponding unknown 
parameters have been fixed from the triton binding energy and nd scattering length, 
which allowed us to make parameter-free predictions for various other observables. For 
the a-particle binding energy we find at NNLO 

BE( 4 He) = -29.51 ... - 29.98 MeV for A = 500 . . . 600 MeV, 

which is rather close to the empirical value. Notice that the cut-off dependence is strongly 
reduced compared to the NLO result. One also observes improvement for various 3N 
scattering observables when going from NLO to NNLO, see Fig. |2]for two examples. 

3. IMPROVING THE CONVERGENCE OF THE CHIRAL EXPANSION 
FOR NUCLEAR FORCES 

As already pointed out in the previous section, the numerically large values of the LECs 
q found in the 7rN system lead to an unphysically strong attraction of subleading TPE. 
Kaiser et al. have performed parameter-free perturbative calculations of the peripheral 
NN partial waves jE] and demonstrated that the nice agreement with the data in D- and 
F-waves at NLO is destroyed after accounting for the NNLO TPE contribution, which 
might indicate problems with the convergence of the chiral expansion. Although these 
problems have been avoided in the previously discussed analysis of the 2N, 3N and 4N 
systems by a reduction of the values of c^a, the situation is clearly not satisfactory, since 
the reduced values of these LECs are not compatible with the ttN system. 

In the following I will explain the origin of the slow convergence and present a new 
method to improve it. To be specific, consider the isoscalar central part of the subleading 
TPE which results from the triangle diagram and is given by 

™ = H / A ^-tr+^Lr?^ ^ + ^ " «"» • (1) 

where q is the nucleon momentum transfer and q = \q\, I = The integral is cubically 
divergent and needs to be regularized. Applying dimensional regularization (DR) one 
finds: 

Mq) = -iH* (2M ' (2ci : cs) ~ csg2) ml + q2) h arctan 2k: + ■ ■ ■ ■ (2) 

2 These values have been corrected to adjust for pp and nn forces missing in the calculations. 
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Figure 3. The potential Vc in r-space. The solid line (shaded band) shows the DR 
(spectral function regularized, A = 500 . . . 800 MeV) result. The dashed (dashed-dotted) 
line refers to the phenomenological a (a + u) + p) contributions based on the isospin triplet 
configuration space version (OBEPR) of the Bonn potential [T7] . 



The ellipses refer to polynomial (in q 2 ) terms of the kind a + f3q 2 . In this section I will be 
interested in D and higher partial waves, where such terms do not contribute. In order 
to obtain the potential in coordinate space one has to make an inverse Fourier-transform 
of Vc{q) in eq. (j2J). The ordinary inverse Fourier-transform is obviously not possible due 
to the the fact that Vc(q) growth with q. One can nevertheless obtain Vc(r) at each r > 
using the spectral function representation jHj: 



V c (q) 



7T 



dfl 



2M-K 



/i 3 p 2 + q 2 



(3) 



where the spectral function p(p) can be obtained from Vc(q) in eq 



via 



3 [V c (0 + - ip)] = -^j {2M 2 (2 Cl - c 3 ) + c 3 /, 2 ) (2M 2 - p 2 )\p - 2M W ) 



(4) 



In eq. the twice subtracted dispersion integral is given which is needed in order to 
account for the large-/! behavior of p(p). 

The inverse Fourier-transform in terms of the spectral function p(p) can easily be 
evaluated via 



Vn(r) 



2n 2 r 



dp:p:e ^ r p(p). 



(5) 



Substituting p{p) from eq. (jlj) into eq. © and using for the LECs Ci j3 the central values 
from ref. [TSJ, c\ = —0.81 GeV -1 and C3 = —4.70 GeV -1 , one obtains the coordinate space 
representation of the potential Vc(r) shown by the solid line in Fig. El The central part 
of the NNLO TPE potential turns out to be several times more attractive at intermedi- 
ate distances than the phenomenological a (a + uj + p) contributions. This unphysical 
attraction shows up in the D- and F-wave phase shifts as depicted in Fig. which have 
been calculated using the Born approximation. While phase shifts at LO (based on OPE) 
and NLO (based on OPE + leading TPE) are in a fair agreement with the data, NNLO 
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Figure 4. Peripheral NN phase shifts. Dotted line refer to LO, dashed and solid lines (light 
and dark shaded bands) to NLO and NNLO calculations using DR (the new regularization 
scheme). Solid dotes are the results from Nijmegen phase shift analysis. 



results (based on OPE + leading TPE + subleading TPE) disagree significantly from the 
data. 

The origin of the unphysical attraction at NNLO can be traced back by looking at 
the integral in eq. (|5|). In Fig. 5 the integrand in eq. (jSJ is plotted versus /a at large 
(r = 2M~ l ), intermediate (r = M~ l ) and short (r = 0.5M" 1 ) distances. While at large 
distances the integral is dominated by 1ow-/a (of the order p ~ 350 MeV) components, 
already at intermediate distances rather high-/! (of the order p ~ 600 MeV) components 
in the spectral function provide a dominant contribution. Clearly, at shorter distances 
even higher-/! components become important. Chiral EFT can hardly provide convergent 
results for the spectral function at p ~ 600 MeV and higher. Instead of keeping such large- 
p components in the regularized loop integral expressions it is advantageous to perform 
the integration in eq. (j3J) only over the low-/x region, where chiral EFT is applicable. This 
can be achieved by introducing the regularized spectral function 

p(p)^p\p)=p(p)9(\-p), (6) 

with the reasonably chosen cut-off A < M p . Certainly, taking a too small A in eq. (JHJ) will 
remove the truly long-distance physics while too large values for the cut-off may affect 
the convergence of the EFT expansion due to inclusion of spurious short-distance physics. 
Notice that a very similar idea with the finite momentum cut-off has been successfully 
applied to improve the convergence of the SU(3) baryon chiral perturbation theory j23 25J. 

In Fig. El we show Vc(r) obtained using the spectral function regularization eq. (JOJ) with 
A = 500 . . . 800 MeV. The strongest effects of the cut-off are observed at intermediate 
and short distances, where the unphysical attraction in dimensionally regularized TPE is 
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Figure 5. The (normalized) 
integrand J(/i) in Eq. (JHJ) for 
different distances r. 
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greatly reduced. On the other hand, the asymptotic behavior of the potential at large r 
is not affected by the choice of regularization. The results for D and F-waves are greatly 
improved when the spectral function regularization is used instead of dimensional one, as 
documented in Fig. El The error of about 10 (1)° at E l&h = 300 MeV for the D (F-) 
waves appears reasonable. 

In ref. |22] we have demonstrated that the spectral function regularization is equivalent 
to (finite) momentum cut-off regularization of pion loop integrals. It should be understood 
that this new regularization scheme does not introduce any model dependence in the EFT 
procedure as soon as A is chosen of the order of (or larger than) M p . Various choices for 
A (including A = oo, which is equivalent to dimensional regularization) differ from each 
other by higher-order contact terms and lead to exactly the same result for observables 
provided one keeps terms in all orders in the EFT expansion. In ref. [22] we have also 
demonstrated how to systematically perform renormalization in this scheme, see also 
|23|24| for the related discussion. 

Last but not least, our choice of regularization in eq. © is by no means unique. Different 
choices lead to equivalent results for the potential up to higher order terms and may be 
used as well. The advantage of the form eq. (jHJ) is that it does not generate spurious long- 
range contributions which are suppressed by inverse powers of A [22]- Simple analytical 
expressions for the regularized loop functions at NNLO are given in ref. (22]- 

4. TWO NUCLEONS WITH THE NEW REGULARIZATION SCHEME 

Application of the new regularization scheme to the low partial waves in the non- 
perturbative regime is the last topic I would like to address. The explicit expressions 
for the potential up to NNLO derived with the new regularization scheme are given in 
ref. j2E] • I n that work the central values from the Q 3 -analysis of the ttN system ^H] have 
been adopted for the LECs Ci )4 : C\ = —0.81 GeV -1 , c 4 = 3.40 GeV -1 . For the constant c 3 
the value C3 = —3.40 GeV -1 has been used which is on the lower side but still consistent 
with the results from ref. (C3 = —4.69 ± 1.34 GeV -1 ). Interestingly, similar values for 
this LEC have been extracted recently from matching the chiral expansion of the nucleon 
mass to lattice gauge theory results at pion masses between 500 and 800 MeV [27] . 

Using the effective potential the bound and scattering states are generated by solving 
the Lippmann-Schwinger (LS) equation which has to be regularized because of the in- 
correct behavior of the potential at large momenta (even after removing the large-mass 
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Figure 6. Selected low NN partial waves calculated with the new regularization scheme. 
For notations see Fig. 0J 



components in the spectrum). We regularize the LS equation in a standard way using 
the exponential regulator function f A (p) = exp[— p 6 /A 6 ], see for more details. In that 
reference A has been varied in the range 450. . .600 MeV at NLO and 450. . .650 MeV at 
NNLO. The cut-off in the spectral function has been varied independently in the range 
500. . .700 MeV. 

For any choice of A and A (denoted in ref. by A) the 9 LECs accompanying the 
contact interactions have been fixed from a fit to np low partial waves at low energy 
(-Eiab < 100 MeV) in the same way as in |8l9j . This allows to make predictions for 
higher energies/partial waves. Results for selected partial waves are shown in Fig. [HI The 
description of the data improves at NNLO compared to NLO in nearly all channels, which 
is entirely due to inclusion of the subleading TPE contribution. The improvement is also 
observed for the scattering lengths and effective range coefficients in both S-waves as well 
as for the deuteron properties [26 . Notice that contrary to the results of ref. jH], where 
DR has been applied to calculate pion loop integrals, there are now no spurious deeply 
bound states. 

5. SUMMARY AND OUTLOOK 

I have discussed several topics related to the EFT description of few-nucleon systems. 
At NNLO, the complete analysis of 2N, 3N and 4N systems has been performed including 
for the first time the chiral 3N force. Promising results have been obtained at NNLO for 
various Nd scattering observables as well as for the a-particle binding energy. 

The problem with the unphysically strong attraction of the chiral TPE has been resolved 
using the new regularization scheme based on the spectral function representation. The 
method has already been successfully applied to the 2N system at NNLO, where it allows 
to significantly improve the convergence of the EFT expansion. It will be interesting to 
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reconsider the 3N and 4N systems using this formalism. 

Finally, it is now of utmost importance to investigate the next order (NNNLO) to be able 
to make conclusions about the convergence of the chiral EFT for nuclear systems. Work 
along these lines is in preparation [2E] (for a first attempt using dimensional regularization 



see |22|). 
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